Phase separation versus supersolid behavior in frustrated antiferromagnets 
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We investigate the competition between spin-supersolidity and phase separation in a frustrated 
spin-lialf model of weakly coupled dimers. We start by considering systems of hard-core bosons 
on the square lattice, onto which the low-energy physics of the herein investigated spin model can 
be mapped, and devise a criterion for gauging the interplay between supersolid order and domain 
wall formation based on strong coupling arguments. Effective bosonic models for the spin model 
are derived via the contractor renormalization (CORE) algorithm and we propose to combine a 
self-consistent cluster mean-field solution with our criterion for the occurrence of phase separation 
to derive the phase diagram as a function of frustration and magnetic field. In the limit of strong 
frustration, the model is shown to be unstable toward phase separation, in contradiction with 
recently published results. However, a region of stable supersolidity is identified for intermediate 
frustration, in a parameter range not investigated so far and of possible experimental relevance. 

PACS numbers: 75.10.Jm, 03.75.Nt, 05.30. Jp 



I. INTRODUCTION 

Dimer-based antiferromagnets (DAFs) under a mag- 
netic field are promising candidates for displaying new 
phases of bosonic matter.^ Magnetic excitations in such 
systems, termed triplons, are well described by lattice 
models of interacting bosons, whose density can be finely 
tuned by varying the magnitude of the applied field. ^''^ 
Experimentally, field-induced Bose-Einstein condensa- 
tion (BEC) of triplons has been observed in a number 
of DAFs (see the review Ref. 1) and, remarkably, ex- 
otic quantum criticality has been detected in the spin- 
dimer compound BaCuSi206.'*"^ The presence of mag- 
netic frustration further adds to the rich phenomenol- 
ogy of these systems by enhancing repulsive interac- 
tions between triplons, something that may eventually 
stabilize incompressible phases that break the lattice's 
translational symmetry.^'* Such crystalline phases are 
for instance realized in the Shastry-Sutherland material 
SrCu2(B03)2,^'^° where they are signaled by a series of 
magnetization plateaux at unconventional fillings stabi- 
lized by complex triplon interactions.** 

The occurrence of both BEC and solid phases in the 
phase diagram of DAFs under magnetic field suggests 
that the magnetic equivalent of the phase simultaneously 
displaying diagonal and off-diagonal order known as su- 
persolid (SS)*^"** may be realized in these systems. In- 
deed, insofar as more exotic possibilities are excluded,*^ 
according to the Ginzburg-Landau- Wilson paradigm a 
continuous transition between phases breaking different 
symmetries (as it is the case with BEC and crystalline 
phases) is precluded and we are therefore left with two 
possibilities: (i) a first-order transition or (ii) the ap- 
pearance of an intermediate phase, where both order 
parameters coexist, termed spin-supersolid (spin-SS) in 



the present context. The latter possibility has been 
first investigated by Momoi and Totsuka for the Shastry- 
Sutherland model in the vicinity of half- and third-filling 
plateaux,*^ based on a mean- field analysis of an effec- 
tive bosonic model derived up to third-order in the inter- 
dimer coupling. More recently,** state-of-the-art tech- 
niques have been employed in deriving effective models 
that improve upon the third-order effective Hamiltonian 
of Ref. 16. Unfortunately, the reliability of these meth- 
ods is still limited to inter-dimer couplings equal to, at 
most, one-half of the intra-dimer coupling, and in that 
parameter range the different plateaux seem to be sep- 
arated by first-order transitions without any convincing 
evidence of spin-SS phases.** 

The situation is much clearer for the DAF investi- 
gated in Ref. 17 where repulsion among triplons are en- 
hanced due to the strong Ising-like character of the inter- 
dimer exchange [see Eq. (1) in Ref. 18], making room for 
checkerboard solid (CBS) and spin-SS phases to emerge. 
The absence of frustration allows for quantum Monte 
Carlo (QMC) simulations to be performed and, in this 
way, the occurrence of a spin-SS phase for the model 
studied in Ref. 17 has been firmly established. However, 
such a strongly anisotropic Hamiltonian is unrealistic for 
Mott insulating materials and further investigations of 
models where the kinetic energy is instead reduced by 
frustration of isotropic couplings*^'^° are clearly called 
for if connection to experiments is ever to be made. 

In this context, the recent report of a spin-SS phase in 
a spin-half frustrated DAF by Chen et al.,'^'^ who have 
relied on a novel tensor-product algorithm, is an impor- 
tant result. However, in view of the first-order transitions 
observed in the related case of the Shastry-Sutherland 
model, a systematic investigation of the possibility of 
phase separation (PS) is still required. 
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FIG. 1: (Color online) (a) A DW between mismatching do- 
mains in a CBS doped with holes is highlighted: open cir- 
cles represent hard-core bosons (triplets) and holes/singlets 
in the upper (lower) domain are indicated by upward (down- 
ward) triangles; doped holes are shown as light-filled upward 
or downward triangles, (b) A simplified model for the DW, 
valid for Vi/ti ^ 1, is defined on a "comb" geometry: holes 
hop (with amplitude t\) through the links indicated by dashed 
lines and repel, with strength Vi, one another along the ver- 
tical nearest-neighbor links indicated by solid lines. 



1.55 

1.5 
1.45 

1.4 
1.35 

1.3 
1.25 

1.2' 
1.15 



10 



A /I , = 3 

h 

O n=6 



III 
10-^ 




10^ 



t./V , 



III 
10" 



10 



10 



FIG. 2: (Color online) c = —Eow/nhti obtained from EDs 
on the geometry depicted in Fig. 1(b), for clusters comprising 
from iV = 8 (ni, = 2) to iV = 24 (rih = 6) sites. 



In this paper, wc investigate the interplay between SS 
order and PS in the frustrated DAF analyzed in Ref. 20. 
We begin by estimating the energetic gains behind PS 
and supersolidity for hard-core bosons on the square lat- 
tice by relying on strong coupling arguments, and intro- 
duce an indicator of the instability toward PS. We then 
proceed to the analysis of effective bosonic models ob- 
tained from the application of the Contractor Renormal- 
ization (CORE) algorithm^^ to the DAF investigated in 
Ref. 20. The so-obtained effective Hamiltonians are stud- 
ied by performing self-consistent cluster mean-field the- 
ory (SCMFT) calculations and tendency toward PS is 
gauged through means of the aforementioned indicator. 



II. PHASE SEPARATION VERSUS 
SUPERSOLIDITY FOR LATTICE BOSONS 

In this section, we analyze the interplay between PS 
and supersolidity in models of hard-core bosons on the 
square lattice, onto which the low-energy physics of the 
spin model considered in the remainder of this paper can 
be mapped. 



A. Instability to Domain Wall Formation 

We start by considering the simplest model of hard- 
core bosons on the square lattice, the so-called t — V 
model: 
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(1) 

n.i — h\b.i is the occupation number operator for holes^^ at 
site i (a hard-core constraint is imposed), whose density 
is controlled by the chemical potential fi. Here, only 



nearest-neighbor (NN) hopping (ti) and interaction (Vi) 
terms arc considered but, as will be discussed later, the 
effective model derived for the spin model Eq. (3) further 
comprises longer-ranged and/or multi-body couplings. 

The model Eq. (1) had been for many years conjec- 
tured to support SS phases (see Refs. 23,24 and references 
therein), but a more systematic numerical analysis^^ later 
disproved earlier evidences in favor of this scenario and 
showed that, instead, PS prevents the occurrence of su- 
persolidity. An intuitive explanation for such behavior 
was put forward in Ref. 25 by relying on strong coupling 
arguments. Following this analysis, ^^ holes [or also parti- 
cles in the case of Eq. (1);^^ however, the effective models 
to be analyzed later lack particle-hole symmetry] doped 
into the CBS ground-state of Eq. (1) for Vi/ti » 1 at 
half-filling would delocalize with an effective hopping am- 
plitude proportional to tf/Vi and eventually condense, 
giving origin to SS order. However, this last conclusion 
is fiawed in that it ignores the possibility of PS. Indeed, 
in the strongly interacting regime domain wall [DW, de- 
picted as the shaded region in Fig. 1(a)] formation is 
energetically favored for the model Eq. (1): the ener- 
getic gain per doped hole (we denote the number of doped 
holes by ?7h) is linear in ti under these circumstances, 
E-Dw/rih ^ ~cti with c G [1,2].25 

Since one of our primary goals in the present work is 
to investigate the interplay between PS and supersolidity, 
so to be able to decide which among the two possibilities 
take place for the spin model herein analyzed [Eq. (3)], 
we would like to obtain a more accurate estimate for 
i^Dw/'^h^i; in other words, we would like to pinpoint 
the actual values assumed by c £ [1,2]. In achieving this 
goal we completely ignore fluctuations in CBS-ordered 
regions away from DWs, a supposedly good approxima- 
tion for Vi/ti ^ 1, and consider a simplified t — V-like 
model at half-filling (quarter-filling for doped holes) de- 
fined on the "comb" geometry depicted in Fig. 1(b). In 
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FIG. 3: (Color online) (a-b) Correlated hoppings behind the 
"leapfrog mechanism" for supersolidity [holes hop in between 
red and light-blue sites only if the dark circles are occupied 
by holes — in (a), at least one of the sites must be occupied; 
if both are, the amplitude is 2si], that allow extra holes to 
delocalize in a CBS background by leapfrogging on the other 
sublattice (c). Adapted from Ref. 18. 



such a simplified model hopping processes with ampli- 
tude f 1 only take place in between NN sites linked by the 
com.b's "teeth" and interaction Vi is only active for holes 
sitting on NN sites along the "backbone" [see Fig. 1(b)]. 
In Fig. 2 we plot results for c = —Euw/nhti, as a 
function of ti/Vi, obtained from exact diagonalizations 
(EDs) of the just discussed simplified model on the comb 
geometry depicted in Fig. 1(b), for clusters comprising 
up to iV = 24 sites (thus, up to rih = 6 doped holes). We 
first notice that c w 1.2 in the limit of large Vi/ti. On 
the other hand, since fluctuations away from the DW are 
ignored in our analysis, we expect our ED results to un- 
derestimate c for small values of Vi/ti. Nonetheless, we 
keep this limitation in mind and throughout the remain- 
der of this paper rely on ED results in estimating the DW 
energy even in the weakly interacting regime Vi/ti '^ 2. 
In doing so, we take advantage of the very small finite 
size effects in the data shown in Fig. 2 and obtain c from 
EDs on finite clusters. 



B. Leapfrog Mechanism for Supersolidity 

In a previous work,^® we have shown that the insta- 
bility toward PS is suppressed in models of hard-core 
bosons on the square lattice that include, in addition to 
the terms comprised in the t — V model [Eq. (1)], the 
correlated hopping processes with amplitudes si and S2 
depicted in Fig. 3. Indeed, such correlated hopping terms 
have been shown^^'^^ to favor supersolidity by allowing 
doped holes to delocalize on top of a CBS by "leapfrog- 
ging" on the other sub-lattice. We wish now to devise 
a criterion for determining how large should the ampli- 
tudes si and S2 be so to inhibit DW formation and thus 
stabilize a SS. In doing so, we once more rely on strong 
coupling arguments, and estimate the energetic gain as- 
sociated to the leapfrog processes represented in Fig. 3 
as 



^ss/mi = -4(2|.si| + |s2|) 



(2) 



That is, Egg/nii is simply the ground-state energy of a 
single hole doped into a "frozen" CBS, an approximation 




FIG. 4: (Color online) (a) Antiferromagnetic bilayer investi- 
gated in this paper [Eq. (3)], with couplings: J± (thick ver- 
tical lines), J|| (thiner in-layer lines) and Jx (dashed lines). 
A magnetic field h promotes singlets (vertical pairs of filled 
circles) to triplets (pairs of open circles); a CBS configuration 
at half-filling is depicted, (h) N = 2 x 2 cluster for SCMFT: 
interactions (thick black lines) involving only in-cluster sites 
(dark-filled circles) are treated exactly while couplings to the 
environment (grey lines) in a MP way. Although only NN 
bonds are depicted, the effective model from CORE also in- 
cludes longer-ranged terms. 



expected to hold for Vi/ti ^ 1, hopping via the processes 
with amplitude si and S2 [Fig. 3(a-b)]. As it happens for 
our estimate ExysN/ni^ ~ —cti obtained in Sec. II A, we 
expect Ess/n-h as given by Eq. (2) to underestimate the 
actual energetic gain associated to the leapfrog processes. 

We combine the just presented analysis and the one 
discussed in Sec. II A concerning DW formation and in- 
troduce an indicator for analyzing the interplay between 
supersolidity and PS: the difference between £'ss/"-h 
[Eq. (2)] and i?Dw/"-h = — cti (obtained from EDs by us- 
ing Vi/ti as input; see Sees. II A), our estimates for the 
energetic gains respectively associated to each of these 
possibilities. Since both estimates are obtained from 
strong coupling analysis, the indicator (i?ss ^ -E'Dw)/"h 
can only be expected to be accurate in the limit of 
Vi/ti ^ 1. However, we keep this limitation in mind 
and in the analysis to be performed in Sec. Ill, we rely 
on (-Ess ~ £'Dw)/?^ii as an indicator even for couplings 
Vi/ti - 2. 



III. FRUSTRATED SPIN MODEL 

A. Model and Effective Hamiltonian 

In most of the lattice models that have, so far, been 
shown to display SS properties the effective repulsion 
(necessary to destabilize the uniform supcrfluid and in- 
duce a SS state) stems from the presence of nearest- 
neighbor repulsive terms in XXZ Hamiltonians.^^^^^'^^ 
While these anisotropic models are interesting from a 
theoretical perspective, their strong anisotropic charac- 
ter renders them unrealistic for antiferromagnetic Mott 
insulators. More promising in this sense is the frustrated 
spin-half Hamiltonian analyzed in Ref. 20, defined on a 



bilayer geometry [Fig. 4(a)] 
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(i, j) denotes NN sites in each square layer a of the frus- 
trated bilayer depicted in Fig. 4(a). J_l couples spins in 
different layers to build the basic dimers of the model 
(we set J_L = 1). The applied magnetic field h acts as a 
chemical potential, promoting spin-dimers from a singlet 
(hole) to a triplet (triplon) state. Effective interactions 
appear as the result of in-layer Jh and frustrating Jx 
antiferromagnetic couplings. We remark that the lattice 
depicted in Fig. 4(a) remains invariant if every other spin- 
dimer is rotated by tt and thus Eq. (3) is invariant under 
the transformation J|| ^-> Jx , with the consequence that 
the phase diagram is symmetric about the line Jh = Jx- 
In studying the model of Eq. (3), we adopt an ap- 
proach similar to the one employed in our previous work 
Ref. 18, where a related unfrustrated model was investi- 
gated and to which the reader is referred for details, ^^ 
and derive an effective bosonic model by relying on 
the CORE algorithm. ^^ We consider spin-dimers con- 
nected by J± as elementary blocks and select the sin- 
glet \s) = -i=[|t|) - |;t)] and the S^ = +1 triplet 

|i+) = Itt) as the block states in the CORE expansion: 
for all parameters in Eq. (3) considered in the present 
work, J|| , Jx G [0, 0.5], this choice is justified by the large 
reduced density-matrix weights associated to such block 
states and by the rapid convergence of effective couplings 
for increasing range in the expansion. ^^ Effective cou- 
plings are derived by diagonalizing clusters of coupled 
dimers and by projecting a matching number of low-lying 
cluster eigenstates onto the basis formed by tensor prod- 
ucts of the retained block states, \s) and |i+}.^^ The ef- 
fective bosonic Hamiltonian thus obtained is essentially 
identical to the one derived for the anisotropic model 
studied in Ref. 18, only the magnitudes for each coupling 
being different. Similarly, the effective model obtained 
here is not invariant under particle-hole transformation 
and, in particular, amplitudes for "leapfrog processes" 
are non-zero only when holes are involved. From this last 
observation we expect that only hole-doped SS phases 
can be stabilized in the spin-dimer model Eq. (3) and 
conclude that the effective Hamiltonian is more conve- 
niently expressed in terms of hole operators (n^ = blbi 
is the occupation number for holes at the dimer-lattice 
site Ti). We therefore adopt the same notation as in our 
previous work,^*'^^ to which the reader is referred for a 
complete list of single- and multi-hole interactions and 
hopping processes [see Eqs. (8, B1-B5) in Ref. 18]. We 
also remark that the effective Hamiltonian preserves the 
symmetry of the original model Eq. (3) and remains in- 
variant under Ju O Jx . 




0.4 0.6 0.8 



1.2 1.4 1.6 1.8 
h 



FIG. 5: (Color online) SCMFT results for the condensate 
density po [squares, Eq. (4)], CBS structure factor 5'(7r,7r) 
[circles, Eq. (5)] and magnetization density [triangles, Eq. (6)] 
for the effective CORE Hamiltonian for Eq. (3) with cou- 
plings (J||, Jx) considered in Ref. 20: (a) (0.38,0.15) and (b) 
(0.38,0.21). Successive phases for increasing magnetic field 
h are labeled as: spin-gapped (Mo), condensate (BEC), su- 
persolid (SS), checkerboard solid (CBS) and fully polarized 
(Ml). 



B. Mean-Field Approach 

Effective Hamiltonians resulting from CORE are often 
complex and different strategies may be pursued in trying 
to extract physically sound results from them. One pos- 
sibility is the mean-field (MF) theory of Ref. 18, which 
reproduces semi-quantitatively the results of quantum 
Monte Carlo (QMC) simulations^^ for the anisotropic 
spin-dimer model considered therein. However, MF cal- 
culations are known to overestimate the extent of SS 
phases and it would be desirable to include, at least par- 
tially, effects due to quantum fluctuations. From this per- 
spective, the SCMFT'^'^ that partially takes local quan- 
tum fluctuations into account and has been recently ap- 
plied to the t—V model for hard-core bosons on the trian- 
gular lattice'^^ seems particularly well suited for our pur- 
poses. Indeed, the extent of the SS phase in the ground- 
state phase diagram obtained by applying SCMFT to the 
t ~V model on the triangular latticc"^^ compares consid- 
erably better with results from QMC simulations^^ than 
what is found from a more conventional MF approach. '^^ 

SCMFT is applied by diagonalizing the effective CORE 
Hamiltonian on the N = 2 x 2 cluster depicted in 
Fig. 4(b). In setting the cluster's Hamiltonian, in-cluster 
interactions are treated exactly while couplings to the en- 
vironment in a self-consistent way: for instance, a given 
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FIG. 6: (Color online) Jy = 0.38. (a) Extent of CBS (A/icbs) 
and SS (A/iss) phases [maximum minus minimum value of the 
field h leading to the corresponding phase for given param- 
eters {J,Jx)]- (b) Value of the structure factor [Eq. (5)] at 
the CBS plateau. In (a) and (b), symbols indicate results by 
Chen ei aL^° and lines the here obtained results. 



interaction connecting sites r^ and Yj contributes a term 
proportional to niUj for each in-cluster bond [thick black 
Unes in Fig. 4(b)] and with mean-field terms of the form 
[ni{nj) + {ni)nj] for "bonds" connecting the cluster to 
its environment [grey lines in Fig. 4(b)]. At each step, 
the ground-state for the cluster Hamiltonian is calcu- 
lated and expectation values (n,), (6^) at every site r^ 
computed; these are then used in setting the mean fields 
for the next iteration, until convergence is achieved (see 
Ref. 31 for details). In this way, we compute the conden- 
sate density at the point ko = (tt, tt) 



Po = 



N ^^ 



iko-rj 



(b,) 



(4) 



the CBS structure factor (normalized per site) 



j.i 



and the magnetization along the field direction 



(5) 



(6) 



In Fig. 5 we plot these quantities as a function of 
the magnetic field h for couplings considered in Ref. 20, 
{J\l,J^) == (0.38,0.15) and (0.38,0.21). We first notice 
that the overall agreement between our results and the 
data presented in Ref. 20 is remarkably good.'^'^ For the 
least frustrated case of {J^^, Jx) = (0.38,0.15) [Fig. 5(a)], 
the system first undergoes a quantum transition from 
a spin-gapped (equivalent to a trivial bosonic Mott in- 
sulator with zero-filling for triplons, Mg in Fig. 5) to a 
BEC phase at the lower critical field hd, and then from 



the BEC to a fully polarized phase (trivial Mott insu- 
lator with unitary triplon filling. Mi) at the upper crit- 
ical field hc2- More interestingly, additional CBS and 
SS phases are stabilized for the more frustrated case of 
{J\l,J^) = (0.38,0.21) [Fig. 5(b)]. The existence of a SS 
phase at the low-field boundary of the CBS plateau, with 
finite values for both po ^nd S'(7r,7r) is at least partially 
due to the presence of correlated hoppings for holes in 
the effective CORE Hamiltonian. Indeed, no SS phase 
is observed for an "effective model" obtained by setting 
si = S2 = [Fig. 3(a-b)] while keeping all the other 
effective couplings unchanged. This situation is to be 
contrasted with the first-order transition from CBS to 
BEC at higher fields, explained by the vanishing ampli- 
tudes for correlated hoppings for triplons for all values 

J||,Jx e [0,0.5]. 

We proceed by varying the frustrating coupling Jx 
while fixing Jy — 0.38. In Fig. 6(a) we plot the ex- 
tent of the SS, CBS phase, respectively A/iss,CBS = 

l^max _ l^min Il^max ( umin \ HprintpQ flip imnpr 

"SS.CBS "ss,CBS ["ss,CBS l"ss,CBsJ Qsnoxes tne upper 
(lower) boundary of the SS, CBS phase] as a function 
of Jx/J|[. In order to further gauge the accuracy of the 
here employed CORE-SCMFT approach in Fig. 6(a) our 
results for Ahss and A/icbs are compared against those 
from Ref. 20 and in Fig. 6(b) we plot both our results and 
those from Ref. 20 for the structure factor S{tt, tt) at the 
CBS plateau. Excellent agreement is found in both cases 
and we further remark that our results for S'(7r, tt) in the 
CBS phase in Fig. 6(b) confirm that quantum fiuctua- 
tions are indeed partially taken into account by SCMFT: 
in contrast to what happens with the semi-classical MF 
approach employed in Ref. 18, here the value of ^(Tr, tt) at 
the plateau is somewhat reduced from its classical value 

'5classical(7r,7r) = 1/4. 

At this point, and despite of its aforementioned attrac- 
tive features, it is important to have in mind an impor- 
tant limitation of the here employed SCMFT procedure: 
since calculations rely on diagonalizations of a 2 x 2 clus- 
ter [Fig. 4(b)], only homogeneous solutions, displaying 
order consistent with at most quadrupling of the unit 
cell, are obtainable. This excludes inhomogeneous solu- 
tions such as those associated with the presence of DWs 
[Fig. 1(a)] and has the consequence that our combined 
CORE-SCMFT approach is insensitive to the occurrence 
of PS. In what follows, we rely on the strong coupling 
analysis presented in Sec. II and analyze the interplay 
between PS and spin-SS order in the phase diagram of 
the model Eq. (3). 



C. Phase Diagram and Phase Separation 

We now turn our attention to the obtention of a global 
J|| — Jx phase diagram that may guide the experimental 
search for realizations of spin-supersolidity and therefore 
extend our analysis by varying Jn in Eq. (3). In Fig. 7(a) 
we plot A/iss as a function of Jy, Jx S [0, 0.5], obtained 
from the combined CORE-SCMFT procedure. These 



results suggest that, far from being a rare occurrence, 
spin-supersolidity is widespread throughout the param- 
eter space and can extend over fairly wider ranges of h 
than it is observed for the value Jy =0.38 [couplings con- 
sidered in Ref. 20 are highlighted in Fig. 7(a)]. However, 
under the light of our discussion in Sec. II concerning PS 
in systems of hard-core bosons on a lattice, some caution 
is required in drawing conclusions from the results shown 
in Fig. 7(a). 

Following the discussion in Sec. II, we evaluate {Ess — 
£'Dw)/"-h^i, our indicator for analyzing the interplay be- 
tween PS and supersolidity for hard-core bosons on the 
square lattice, as a function of Jy , Jx £ [0, 0.5] for Eq. (3). 
Ess/n-ti, our strong coupling estimate for the energetic 
gain associated to occurrence of SS order, is readily ob- 
tained by plugging the amplitudes for the leapfrog pro- 
cesses si and S2 [Fig. 3(a-b)] obtained from CORE into 
Eq. (2). On the other hand, the estimate i?Dw/"-h for the 
energy associated to PS is obtained from numerical EDs 
for the simplified model for DWs defined on the "comb 
geometry" discussed in Sec. II A by using the effective ra- 
tio Vi/ti, as obtained from the CORE expansion for each 
set J|[, Jx in Eq. (3), as an input. EDs are performed on 
a small cluster comprising TV = 16 sites (nh = 4 doped 
holes): as mentioned in Sec. II A, this is justified by the 
absence of sizable finite size effects for the data displayed 
in Fig. 2. 

We plot {Ess — EY)w)/n\-iti as a function of Jy, Jx G 
[0,0.5] in Fig. 7(b) and assume that two conditions must 
be simultaneously fulfilled for SS phases to exist for the 
model Eq. (3): (i) a SS must be observed within CORE- 
SCMFT and {li) {Ess - E^^) / n^h < 0. Values of Jy , Jx 
leading to (i?ss — ^Dw)/'^h^i = 0, the threshold value for 
a SS phase to appear, are indicated by the thick continu- 
ous curve in Fig. 7. We notice that not all values of Jy , Jx 
yielding a spin-SS phase within our CORE-SCMFT ap- 
proach fulfill {Ess — £'Dw)/"'hii < and expect PS to 
take place under these circumstances instead. Despite of 
the fact that the condition (-Ess~£'Dw)/^hii < consid- 
erably shrinks the size of the region expected to support 
SS phases from a pure CORE-SCMFT analysis, super- 
solidity is still observed for a wide range of couplings in 
Eq. (3) [Fig. 7(a)], possibly realizable in real magnets. 

Intriguingly, we notice that the parameters (circles in 
Fig. 7) for which a spin-SS phase has been detected 
by Chen et al.,^° and also by the pure CORE-SCMFT 
analysis devised here (Figs. 5 and 6), fail to satisfy 
{Ess - -BDw)/"-h^i < [Fig. 7(b)]. Although we cannot 
exclude the possibility that our criterion, that rigorously 
applies only in the limit Vi/ti 3> 1, is too stringent for 
the frustrated model Eq. (3), we remark that a SS phase 
is obtained within our CORE-SCMFT approach^"' and 
in Ref. 20 even for couplings J^^ ~ Jx , where our strong 
coupling arguments become accurate [the ratio Vi/ti di- 
verges toward the line Ji\ = Jx ; see Fig. 7(b), where con- 
tour levels for Vi/ti are plotted as thin continuous lines]. 
This inconsistency suggests that both CORE-SCMFT 
and the novel algorithm employed in Ref. 20 arc insensi- 




FIG. 7: (Color online) (a) SCMFT results for the extent of the 
SS phase A/iss (see main text) for the frustrated DAP Eq. (3). 
The symmetry Jy -H- Jx has been explored in obtaining the 
data. Regions where supersolidity [PS] is expected, where 
(_Bss — Euw)/n'h_t\ < [(-Ess — Euw)/nt,ti > 0] are marked 
by the label SS [PS]. Dashed lines indicate threshold values 
for a CBS/SS to appear at the mean-field level, (b) (i?ss — 
E-Dw)/nhti, as obtained from EDs on an A = 16 (rih = 4 
doped holes) site cluster with the comb geometry depicted 
in Fig. 1(b), for the model Eq. (3). Contour levels for Vi/fi 
(obtained from the CORE expansion) are indicated by thin 
lines and the values Vi/ti = 2, 4, 6 and 8 are highlighted. In 
both panels, circles indicate couplings investigated by Chen et 
al?^ and the thick line couplings yielding the threshold value 
(£'ss - £'Dw)/nhii = 0. 



five to the instability toward PS in systems of hard-core 
bosons on the square lattice and that the obtcntion of SS 
phase for {Ess — £'Dw)/'^hii > is spurious.^* It would 
therefore be important to further test the ability of the 
algorithm employed in Ref. 20 to detect PS in bosonic 
lattice models by, for instance, checking how it compares 
to QMC for the unfrustrated model studied in Rcfs. 17,18 
regarding this issue. 



IV. SUMMARY 

Summarizing, we have studied a spin-half frustrated 
bilaycr model by combining CORE and SCMFT. Our 
results reveal the presence of a spin-SS phase under ap- 
plied magnetic field, which appears at the edge of a half- 
saturated magnetization plateau and is stabilized by a 
"leapfrog mechanism".^® We address the interplay be- 
tween supersolidity and instability toward PS, that pre- 
cludes the emergence of spin-SS phases, by devising a 
quantitative criterion based on strong coupling argu- 
ments. This criterion is genetically applicable to sys- 
tems of hard-core bosons on the square lattice, and it 
would be interesting to further assess its validity by in- 
vestigating models where the interplay between PS and 
SS can be independently analyzed. By relying on this 



criterion, we obtain a global phase diagram for the frus- 
trated spin-dimer antiferromagnet considered herein, and 
show that a spin-SS phase is stable against PS for cou- 
plings realizable in real magnets. We expect that our 
results may guide the experimental search for systems 
exhibiting spin-supersolidity. 
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